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! Abstract 

Recent developments in heavy quarknoium physics are reviewed in brief, in- 
' eluding (i) nonrelativistic QCD(NRQCD), (ii) the importance of color-octet 

in ' 

(QQ) components in the decay and the production of a physical heavy quarko- 

' nium state, (iii) limitation of the NRQCD factorization and (iv) the double 

OO ' 

<3\ ' humps in the tttt spectrum in the decay T(3S) — * T(lS)mr. 



I. INTRODUCTION 



Heavy quarkonium physics was one of the early applications of perturbative QCD 
(PQCD). In the original treatment of the decay and the production of a heavy quarko- 
^ ■ nium state H, it was assumed that (QQ)i inside if is in a definite 2S+1 Lj color-singlet state. 

And production and decay rates of a heavy quarkonium H are assumed to be factorized into 
(i) the short distance (SD) parts which are calculable in PQCD (in a s (Mg)), and (ii) the 
long distance (LD) parts that may be parametrized in terms of the (QQ)i wave function 
and its derivatives at the origin Jl|. 

This factorization hypothesis works well for the S— wave and P— wave states to lowest 
order in a s . For example, the Xco state decays into light hadrons (LH) through Xco 99 at 
the parton level : 

F(X- - LH) = 1^(0)1". (1.1) 
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However, such a factorization in the CSM fails if PQCD correction is included, which is 
signaled by the infrared divergence in its decay rates ||. This infrared divergence implies 
the failure of the factorization in terms of a single wave function of (QQ)i in H. This 
problem remained unsolved until 1992, when Bodwin et al. pointed out the importance of 
the color-octet (ccjg component of H (see (5)) ||. 

In this decade, the prompt J/if) production at the Tevatron could be measured with the 
development of high resolution vertex detector. The data was astonishing in that the CSM 
underestimates the data by a factor of ~ 10 for J / if) production, even worse for if) production 
(by ~ 30). Possible contributions in the CSM are (i) gluon-gluon fusion (gg — > J/if> + g), 
(ii) gluon fragmentation (g — > J /if) + gg), (iii) c— quark fragmentation (c — > J /if) + c), and 
(iv) cascade decays from P— wave charmonia (x c j(1P) — > J /if) + 7), which is the most 
important for the J / if) productions. However, for the case of if) production at the Tevatron, 
one observed that a exp ~ 30 x cr tll based on CSM. Since there are no known P— wave states 
(such as Xcj(2P)) that can decay into if) , two options are avaiable : (i) hypothetical x c j(2P) 
decaying into if) + 7 with a suitable branching ratio ||], and can be tested in B and T decays 
0, and (ii) color-octet gluon fragmentation : g — > (cc)g followed by (cc)g — ► J /if) + soft gg 
H . This second idea is the main theme of this talk, and the next two sections will be devoted 
to this issue. 

Another less understood phenomenon in the heavy quarkonium physics in 1980's is the 
double humps in tctt spectrum in T(3S) — > T(lS')7r7r []. This issue is related with chiral 
dynamics of pions, rather than heavy quarknoium itself, and will be separately discussed in 
Sec. IV. 

In this talk, we review these recent developments in brief. In Section II, the NRQCD (and 
factorization formula) is introduced. In Section III, we discuss the color-octet mechanism 
for J /if) productions in various high energy processes. In Section IV, the double peaks in 
T(35) — > T(15')7r7r is discussed, and the summary is given in Section V. 

II. NRQCD 

A. NRQCD Lagrangian 

A heavy quarkonium state H is associated with several different scales with the following 
hierarchical structure : 

A Q cd ^ M Q v 2 « M Q v « M Q . (2.1) 

In brief, Mq is the heavy quark mass scale that is a typical energy scale for decay and 
production of H, and PQCD in « s (Mq) becomes applicable. The typical momentum scale 
of Q inside H is p ~ MqV ~ 1/size, and the typical kinetic energy scale of Q inside H is 
K.E. ~ Mqv 2 ~ the level splittings. Finally, the typical energy scale for nonperturbative 
dynamics of light quarks and gluons is characterized by Aqcd- From the observed spectra of 
charmonium and upsilon families, one gets v 2 ~ 0.3 for the if) = (cc) system, and ~ 0.1 for 



although this puzzle has nothing to do with the limitation of the CSM 
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the T = (66) system. For sufficiently heavy Q, one has v ~ a s (Mgt>) > q; s (Mq). Formally, 
the v 2 expansion is more important than a s correction [fjj. 

In order to obtain the NRQCD lagrangian, one integrates out the modes with momentum 
> Mq (no Q/Q creation/annihilation possible), and expand the resulting effective action in 
terms of heavy quark velocity. Thus NRQCD lagrangian posseses SU(3) C gauge symmetry, 
P (parity) and C (charge conjugation) symmetry, rotational symmetry, and approximate 
heavy quark spin symmetry. It is written in terms of heavy quark (-0) and heavy antiquark 
fields (x), and their covariant derivatives : 

£nrqcd = ^ (iD + ^-j ^ + X ] (iD - X + £%ht + 5C (2.2) 

Gluon exchange with k ~ p ~ MqV is included in terms of static potential, A , whereas soft 
gluons with k ~ Mqv 2 are treated in QCD multipole expansion. 5C includes the correction 
terms of 0(Mqv 4 ) and higher. One can determine the coefficients of operators in NRQCD 
lagrangian by calculating the same process in QCD and NRQCD, and expand both results 
in powers of k/MQ ~ v, and match to the desired order in v and a s . 

In NRQCD, a physical heavy quarkonium state is represented as a superposition of 
various (<50)i,8 an d dynamical gluons : \J/if>) = 0(i)|(cc)i[ 3 Si]) + 0(v )\(cc) 8 [ 3 Si]g) + 
0(v 2 )|(cc) 8 [ 1 S'o]) + .... Relative importance of various Fock states in the above equation 
(and NRQCD matrix elements in factorization formulae discussed in the next subsection) 
are determined by velocity scaling laws that can be derived from field equation of motion of 
NRQCD 0. 

B. NRQCD factorization 

Decays of a heavy quarkonium into light hadrons are described in terms of local 4-fermion 
operators in NRQCD lagrangian by optical theorem : 

T{H -> LB) = 21m (B\8jC^ ieTmioa \H) = e J^^ (H\O n (A)\B) (2.3) 

M Q 

This is a double expansion in v 2 (Non Pert.) and a s (Mg) (Pert.). 

One can determine the NRQCD matrix elements either from experimental data |Q or 

lattice QCD 0. For Xco decay, the CMS prediction, (1), is modified (in the NRQCD 
formalism) into 

r( Xc0 - LH) = 7ra 2 (M c ) M c 5 (O^Po)) + N, ™j^ e) WS,)) (2.4) 

where 

(O 1 ) = |(0|x t (-^-^X c0 )| 2 , (2.5) 

(Os> = (Xco\^°T a X ■ xVT>| Xc0 ). (2.6) 

Second term comes from the annihilation of the color-octet (cc)g (S— wave) in the physical 
Xco state. IR div. in the first term at 0(a 3 ) cacelled by the IR div. in the color-octet matrix 
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element in the second term. Similar cancellations were shown to occur also in the P— wave 
charmonium production ||. 

Inclusive cross section for a heavy quarkonium (H) production 

1 d 3 P 

E*to(l + 2 H(P) + X) = E m „ Cm „(0|Ol|0) (2.7) 

or, 

rfa = S n do-(gg[n] + X)(0|O H (n)|0) (2.8) 

C mn takes into account the SD of order 1/M C or less, and therefore are calculable using 
PQCD in a s (M c ). The matrix element (0|O^JO) is the VEV of a four-ferion operators 
of NRQCD, and is proportional to the inclusive transition probability of the perturbative 
QQ[n] state into the physical heavy quarkonium state H. In terms of NRQCD operators, 

OL = X^l^ V H ^JC nX , (2.9) 

where /C's are product of a spin matrix, a color matrix and a polynomial in the covariant 
derivative, D, and the projection operator V is defined as 

V H = E S \H(P = 0, S)(H(P = 0), S\, (2.10) 

where the sum is over soft hadron states S with E < A. 



III. COLOR-OCTET MECHANISM IN J/if> PRODUCTIONS 

A. i\) anomaly at the Tevatron 

As discussed in the introduction, the color-octet gluon fragmentation g — > (ccJgpSi] 
followed by (cc) 8 [ 3 Si] — > J/ip + (soft gluons) may explain the large excess of the ip at the 
Tevatron 0. Qualitatively, one has 

CSM ~ c^t> 3 (SD suppressed, LD enhanced) (3.1) 
COM ~ a s v 7 (SD enhanced, LD suppressed) (3.2) 

Cho and Leibovich included the color-octet 1 Sq and 3 Pj contributions as well as the color- 
octet 3 Si in terms of three NP parameters, and obtained the following constraints on these 
parameters from the Tevatron data [[K]] : 

(0|0 8 J/ YS 3 )|0> = (6.6 ±2.1) x 1(T 3 GeV 3 

WoTePM + (Q|o. w (%)|o) ={22± 5) x 10 _ 2 GeV3 

M 2 C 3 K J 

It is very important to check the idea of the color-octet mechanism and NRQCD factor- 
ization in other processes. In view of this, it is crucial to observe that the color-octet matrix 
elements (0\O^ ( 2S+1 Lj)\0) are universal, i.e. process- independent. Therefore, one can de- 
termine these matrix elements from a (set of) process(es), and then apply to other processes, 
and test the NRQCD factorization A lot of works have been done in this line. To name a 
few, hadroproduction of h^Px) Q, B -> J/ip + X, Z° -> J/i/j + X pf, e+e" -> J ftp + X 



at CLEO ||13|| , to name only a few that were discussed at this workshop. 
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B. B — ► J /if) + X 

The relevant effective Hamiltonian for this decay is 

H eff = ^ V cb V; q [Ci(cc) v „ A (qh) v „ A + C 2 (ct a c) v „ A (qt a b) v „ A ] , (3.3) 

with Ci(rnb) ~ 0.13 and (^(wft) ~ 2.21 in the leading logarithmic approximation. The CSM 
prediction to the lowest order in a s is 

r ( B ^ + 1). = cf (, + 5) r , (3,) 

with f = \V cb \ 2 ^ MlM c (\ - ^fj ■ Using (0\O( /lp ( 3 Si)\0) determined from J/if> -> 

/ + / _ , we get S(S — > J/?/> + X) csm = 0.23% compared with the most recent CLEO data 
B{B — > J/ip+X) = (0.80±0.08)%. Higher order corrections in a s and relativistic corrections 
are not that important ||. Color-octet 3 Si contribution 



T(b - (c C ) 8 [^!] + s^J/i, + X)= 1 ^ ^ ' Cl 1 + ] r (3.5) 



is enhanced because of large Wilson coefficients f|, C 2 ~ 17Ci. Thus, we get moderate 
increase in the branching ratio to 0.58%. Color-octet X S§ and 3 Pj also contribute at the same 
order of v 2 expansion , dependent on (0|Os / ' / '( 1 5 , )|0) and (0|O^ ( 3 P/)|0). J/^ polarization 
in B — > J/ip + X |TJ]], which may provide another useful constraint on the color-octet matrix 
elements. Similar analyses have been done for B — > \c + X, and i? decays into D— wave 
charmonium state |T3[. Here the avaiable phase space is rather small so that the parton 



model description may be a poor approximation. 



C. J/tp Photoproduction 

The J / ip photoproduction (7 + p — > J/ip + X) is described as a parton level subprocess 
■y + g^J/ip + gm the PQCD and CSM |R|. This process is advocated as a nice probe of 
gluon distribution function inside proton. 

PQCD corrections (in the CSM) to the leading order results have been done by M. 
Kramer [|17]] . The impacts of this radiative correction is that the scale dependence of a s and 
structure functions is reduced. However, PQCD is out of control for z > 0.8 and P% < 1 
GeV 2 at HERA energy [y/s lp = 100 GeV. Therefore one needs cuts, z < 0.8 and P% > 1 
GeV 2 , when one employes the PQCD correction to the J/ip photoproduction. 

Color-octet contributions were considered by various groups [[18] It was shown 

that dT/dz spectrum in the high z > 0.9 region blows up, where z is defined as Ej/^/E y 
in the proton rest frame. This phenomenon was originally taken to be a signal that the 
determination of color-octet matrix elements from the Tevatron data on J/ip productions 
are inconsist with the J/ip photoproduction. 

However, this may not be the case because of breakdown of NRQCD near z = 1 as 
discussed by Beneke et al. [Q. In the lowest order in v 2 , one ignores AM = Mh — 2Mq, 
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i.e. set Mjj^ « 2M C . This is not valid near the phase space boundary, since one probes 
the dynamics in detail : AE < AM, especially, when the matrix element has mainly a 
support near the phase space boundary This is the case for J/ip photoproduction at large 
z(> 0.8 ~ 0.9). In such cases, predictions become sensitive to the AM, or momentum carried 
away by light hadrons during the hadronization. One can summarize this effect in terms of a 
universal shape function, as b — > S7 in HQET. Shape function shifts the unphysical partonic 
boundary of phase space to the hadronic one that is physically more sensible. Therefore, 
the problem with J/ip photoproduction at z > 0.9 is probably less serious. 



IV. vrvr SPECTRUM IN X" Ttttt 



The M n7r spectra in ip' — > J/ipTnr and T' — > Ttttt can be understood in terms of QCD 
multipole expansion and the low energy theorem for pions, which dictates the following 
amplitude : 



M = Ae-e 



q 2 + CMl 



(4.1) 



with q 2 = (P1+P2) 2 = M 2 n . This amplitude predicts a peak at high M nn region in agreement 
with the data from ARGUS and CLEO. However, the double humps in Y(3S)(k,e) — > 
T(lS)(k , e )vr(p 1 )7r(p 2 ) could not be understood with the above amplitude. There are several 
proposals to this phenomenon, but none of them were successful. In Refs. 
amplitude was modified as 



21 , the above 



M = Ae-e 



q 2 + BE%E 2 + CM; 



(4.2) 



and added the phase shift informations for J = 0, S and D waves. Then, the authors 
of Refs. |2T| could fit the double humps, and predict various angular distributions. The 
predictions agree with the newest CLEO data except the COS6 1 * distribution. It turns out 
that inclusion of terms higher in pion momenta, such as (Ei + E 2 )q 2 and so on, improves 



the agreement with the data. More systematic study in ChPT is in progress [22 



V. CONCLUSION 



NRQCD provides a theoretical framework in which one can study the perturbative (in 
a s (Mq)) and nonperturbative (in v 2 ) aspects of heavy quarkonium physics. In particular, 
the role of the \{QQ)&g) Fock state in the heavy quarkonium production and its decay can 
be rigorously formulated in the NRQCD. One can fit the J/ip and ip production rates at 
the Tevatron via the color-octet mechanism. In order to check this idea at other processes, 
PQCD corrections are essential. This part has been calculated only recently by Petrelli et 
al. |[23| . The complete phenomenological analysis including this new result has not been 
done yet, however. Finally, Double hump in M^ spectrum can be explained in terms of 
an amplitude that satisfies the low energy theorem for pions, once the D— wave dipion 
amplitude is properly included. 
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